Diverse computational tools towards the understanding of HIV targets and design of potential drug candidates.

Ph. D. University of KwaZulu-Natal, Durban 2014.HIV/AIDS still remains to be a challenging epidemic infecting millions of individuals worldwide. The morbidity and mortality rates of HIV-infected patients has been well documented over the years. Despite on-going HIV/AIDS research and access to antiretroviral therapy, to date still no cure exists for this deliberating disease. In recent years, computational approaches have emerged as close counterparts to experiments in modern drug discovery process and in understanding complex biological phenomena. An array of in-silico computational techniques were implemented ranging from molecular dynamic (MD) simulations, de-novo design, hybrid structure-based and pharmacophore-based virtual screening, quantitative structure-activity relationship (QSAR), homology modeling, principle component analysis (PCA), residue interaction network analysis (RIN), substrate envelope analysis (SEA), to molecular mechanics and quantum mechanics. The first report (Chapter 4), demonstrated a unique strategy for developing dual acting inhibitors against HIV-1 protease (PR) and reverse transcriptase (RT). The designed targets exhibited binding affinities and dual inhibiting activity comparable to, and in some cases better than, known active reference drugs. The second study (Chapter 5), reported the activity of flexible hydroquinone-based compounds as non-nucleoside reverse transcriptase inhibitors (NNRTIs), as proposed by Bruccoleri, where no experimental or computational work supported his proposal. Results concluded that the novel flexible hydroquinone-based compounds showed improved binding affinity as compared to FDA-approved prototype drugs and more specifically potent potential mutant-resistant NNRT inhibitor activity. The third report (Chapter 6), explored the activity of novel CCR5 antagonists as potential HIV- 1 entry inhibitors. Ten scaffolds were identified as novel CCR5 antagonists or potential HIV-1 entry inhibitors. Furthermore, from the generated atom-based 3D-QSAR model, all of the parameters showed certain reliability and feasible predictability to help us design new and high selectivity CCR5 inhibitors. The fourth study (Chapter 7), explored the atomistic basis of why the M184I single mutation renders complete resistance of HIV-1 RT to lamivudine. Multiple molecular dynamics simulations, binding free energy calculations, principle component analysis (PCA) and residue interaction network (RIN) analyses adequately clarified the effect of the M184I mutation on drug resistance to lamvudine. Results presented in this study verified that M184I mutation decreased drug binding affinity, distorted ligand optimum orientation in RT active site and affected the overall protein conformational landscape. The results also provided some potential clues for further design of novel inhibitors that are less susceptible to drug resistance. In the fifth study (Chapter 8), we identified potential HIV-Nef inhibitors by exploiting the structural features of B9 using an integrated computational tools framework. The top identified hit compounds demonstrated comparatively better binding affinities and relatable binding modes compared to the prototype antagonist, B9. Top identified hits were proposed as new potential novel leads targeting HIV-Nef with a detailed analysis of their respective binding modes. The sixth report (Chapter 9), aimed to reveal the dimer packing and unpacking phenomena of HIV-Nef in its apo and inhibitor bound conformations using molecular dynamic simulations. Results verified a more conformational flexible nature of HIV-Nef dimer in the absence of an inhibitor.as compared to B9 bound conformation of HIV-Nef, which was found to be more conformationally rigid with a lesser inter-dimeric association. We believe that the results obtained from these several studies could be of great benefit in the development of more effective therapeutic interventions for the treatment and cure of HIV/AIDS

Integrated Computational Tools for Identification of CCR5 Antagonists as Potential HIV-1 Entry Inhibitors: Homology Modeling, Virtual Screening, Molecular Dynamics Simulations and 3D QSAR Analysis

Using integrated in-silico computational techniques, including homology modeling, structure-based and pharmacophore-based virtual screening, molecular dynamic simulations, per-residue energy decomposition analysis and atom-based 3D-QSAR analysis, we proposed ten novel compounds as potential CCR5-dependent HIV-1 entry inhibitors. Via validated docking calculations, binding free energies revealed that novel leads demonstrated better binding affinities with CCR5 compared to maraviroc, an FDA-approved HIV-1 entry inhibitor and in clinical use. Per-residue interaction energy decomposition analysis on the averaged MD structure showed that hydrophobic active residues Trp86, Tyr89 and Tyr108 contributed the most to inhibitor binding. The validated 3D-QSAR model showed a high cross-validated rcv2 value of 0.84 using three principal components and non-cross-validated r2 value of 0.941. It was also revealed that almost all compounds in the test set and training set yielded a good predicted value. Information gained from this study could shed light on the activity of a new series of lead compounds as potential HIV entry inhibitors and serve as a powerful tool in the drug design and development machinery

Individual and Population Level Impact of Key HIV Risk Factors on HIV Incidence Rates in Durban, South Africa.

We aimed to estimate the individual and joint impact of age, marital status and diagnosis with sexually transmitted infections (STIs) on HIV acquisition among young women at a population level in Durban, KwaZulu-Natal, South Africa. A total of 3,978 HIV seronegative women were recruited for four biomedical intervention trials from 2002-2009. Point and interval estimates of partial population attributable risk (PAR) were used to quantify the proportion of HIV seroconversions which can be prevented if a combination of risk factors is eliminated from a target population. More than 70% of the observed HIV acquisitions were collectively attributed to the three risk factors: younger age (<25 years old), unmarried and not cohabiting with a stable/regular partner and diagnosis with STIs. Addressing these risks requires targeted structural, behavioural, biomedical and cultural interventions in order to impact on unacceptably high HIV incidence rates among young women and the population as a whole

Integrated Computational Tools for Identification of CCR5 Antagonists as Potential HIV-1 Entry Inhibitors: Homology Modeling, Virtual Screening, Molecular Dynamics Simulations and 3D QSAR Analysis

Using integrated in-silico computational techniques, including homology modeling, structure-based and pharmacophore-based virtual screening, molecular dynamic simulations, per-residue energy decomposition analysis and atom-based 3D-QSAR analysis, we proposed ten novel compounds as potential CCR5-dependent HIV-1 entry inhibitors. Via validated docking calculations, binding free energies revealed that novel leads demonstrated better binding affinities with CCR5 compared to maraviroc, an FDA-approved HIV-1 entry inhibitor and in clinical use. Per-residue interaction energy decomposition analysis on the averaged MD structure showed that hydrophobic active residues Trp86, Tyr89 and Tyr108 contributed the most to inhibitor binding. The validated 3D-QSAR model showed a high cross-validated rcv2 value of 0.84 using three principal components and non-cross-validated r2 value of 0.941. It was also revealed that almost all compounds in the test set and training set yielded a good predicted value. Information gained from this study could shed light on the activity of a new series of lead compounds as potential HIV entry inhibitors and serve as a powerful tool in the drug design and development machinery

Single Active Site Mutation Causes Serious Resistance of HIV Reverse Transcriptase to Lamivudine : Insight from Multiple Molecular Dynamics Simulations

Molecular dynamics simulations, binding free energy calculations, principle component analysis (PCA), and residue interaction network analysis were employed in order to investigate the molecular mechanism of M184I single mutation which played pivotal role in making the HIV-1 reverse transcriptase (RT) totally resistant to lamivudine. Results showed that single mutations at residue 184 of RT caused (1) distortion of the orientation of lamivudine in the active site due to the steric conflict between the oxathiolane ring of lamivudine and the side chain of beta-branched amino acids Ile at position 184 which, in turn, perturbs inhibitor binding, (2) decrease in the binding affinity by (~8 kcal/mol) when compared to the wild-type, (3) variation in the overall enzyme motion as evident from the PCA for both systems, and (4) distortion of the hydrogen bonding network and atomic interactions with the inhibitor. The comprehensive analysis presented in this report can provide useful information for understanding the drug resistance mechanism against lamivudine. The results can also provide some potential clues for further design of novel inhibitors that are less susceptible to drug resistance

Dynamic features of apo and bound HIV-Nef protein reveal the anti-HIV dimerization inhibition mechanism

<div><p></p><p>The first account on the dynamic features of Nef or negative factor, a small myristoylated protein located in the cytoplasm believes to increase HIV-1 viral titer level, is reported herein. Due to its major role in HIV-1 pathogenicity, Nef protein is considered an emerging target in anti-HIV drug design and discovery process. In this study, comparative long-range all-atom molecular dynamics simulations were employed for apo and bound protein to unveil molecular mechanism of HIV-Nef dimerization and inhibition. Results clearly revealed that B9, a newly discovered Nef inhibitor, binds at the dimeric interface of Nef protein and caused significant separation between orthogonally opposed residues, namely Asp108, Leu112 and Gln104. Large differences in magnitudes were observed in the radius of gyration (∼1.5 Å), per-residue fluctuation (∼2 Å), C-alpha deviations (∼2 Å) which confirm a comparatively more flexible nature of apo conformation due to rapid dimeric association. Compared to the bound conformer, a more globally correlated motion in case of apo structure of HIV-Nef confirms the process of dimeric association. This clearly highlights the process of inhibition as a result of ligand binding. The difference in principal component analysis (PCA) scatter plot and per-residue mobility plot across first two normal modes further justifies the same findings. The in-depth dynamic analyses of Nef protein presented in this report would serve crucial in understanding its function and inhibition mechanisms. Information on inhibitor binding mode would also assist in designing of potential inhibitors against this important HIV target.</p></div

Crude incidence rates for HIV seroconversion.

<p>Crude incidence rates for HIV seroconversion.</p